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Abstract 

The dynamics and the transition to the centrifugal regime are studied analytically 
and numerically for particles in rotating drum. The importance of the particle-wall 
friction coefficient is demonstrated by studying first the motion of one non-rotating 
particle where three different regimes are found in the transition to the centrifugal 
motion. When a few rotating particles are considered, they behave similarly to one 
non- rotating particle in the low friction limit. A critical particle number is necessary 
to reach the centrifugal regime for which an analytic expression is derived in the 
limit of negligible inter-particle friction. 

PACS: 81.05.Rm, 46.30.Pa, 05.70.Fh, (46.10+z) 



1 Introduction 

Granular materials can show fascinating behaviour in many experimental situations 0|. 
When placed in a slowly rotating drum, most particles undergo a solid body rotation 
by following the external motion and a thin layer of intermittently flowing particles 
(avalanches) is forming close to the surface. With increasing rotation speed, the separa- 
tion time of avalanches will decrease until no individual avalanches are detectable and a 
nearly constant particle flow is found along the free surface j|] . For even higher rotation 
speeds, the free surface does not stay flat but rather deforms with increasing rotation 
speed. These deformations mostly start from the lower boundary inwards ||. When the 
rotation speed is further increased, the particle motion is dominated by the centrifugal 
force leading to a ring formation close to the outer wall boundary [Q]. In order to under- 
stand the different flowing regimes, it is instructive to study the particle dynamics on a 
microscopic level, i.e. on a particle basis. This is especially helpful for studying different 
frictional forces S which play an important role in many other technological problems 
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as well 0. In a simplified view, one can look at the particles undergoing a solid body 



rotation as sticking to their neighbours whereas the particles in the fluidized surface 
layer undergo mostly sliding collisions. A similar particle behaviour can be observed for 
a particle placed on a moving plate and attached to one side wall via a spring. For low 
surface velocities, a stick-slip motion is observed and above a critical velocity, the stick- 
slip motion disappears and only steady sliding is observed j?], ||. But this setup might 
not allow for long enough observation times in order to investigate transient behaviour. 
We propose another setup, namely to place the particle(s) in a rotating drum. 

The paper is organized in the following way: in section ||, we will sketch the system 
in mind, explain briefly an experimental realization and derive the equation of motion 
for one non-rotating particle in the full-sliding limit. In section ^, we will present a 
numerical model which is not restricted to the full-sliding limit but is capable to explore 
the full particle dynamics including the transition from sticking to sliding motion. This 
will be compared in section Q to the analytic solution and the transition to the centrifugal 
regime will be discussed. In section [EL the collective behaviour of a few rotating spheres 
in the drum will be investigated and a connection to the one particle case can be made. 
The conclusions will round off this paper. 

2 Setup and Theory for One Non-Rotating Particle 

Different friction regimes for particle-wall contacts can be studied by looking at the 
dynamics of a non-rotating particle placed in a rotating drum. The setup is sketched in 
Fig. 0a where the particle is approximated by a sphere shown in grey. The drum rotates 
with a constant angular velocity Q, gravity is pointing downwards and the angle ip of 
the particle's center of mass is measured from the vertical. The distance between the 
drum's and the particle's center is R. In the corresponding experimental system, the 
non-rotating particle was realized by glueing two spheres together and the drum was 
slightly tilted in order to force the particle to slide along one of the end caps || . 

The equation of motion is one-dimensional as long as the particle stays in contact 
with the wall of the drum, i.e. R = constant. For full sliding, \tp\ < f2, the equation of 
motion reads 



where (fo = arctan(// l „) and the Coulomb friction law was used. Here F n denotes the 
normal force and F s the shear force. Having the many particle program at hand, we 
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used MD simulations and compared the results for one non-rotating particle with the 
experiment and the theory. The angle cp as function of time for a typical run is shown 
in Fig. using a drum rotation speed of = 1 Hz and the initial conditions ¥>(0) = 
0(0) = 0. 

Equation can be simplified by the following steps [ 10 1 : setting tp = tp — ipa one 
gets 



Rip = -g Vl+74 sin(v?) + /j, w R<p . (2) 

This equation does not exhibit an explicit time dependence and it is favourable to use 
the transformation u{ip) = <p =^ tp = u'ip = u'u leading to the equation 

Ru'u = -g ^fl + nlj sin(^) + n w Ru 2 . (3) 

By using the transformation v = \v? the equation can be integrated and one obtains 



e -M v (Jp) = V _ rw / e - W gin ^ d ^ + e « v (_^ ) . (4) 

R J- vo 

Please note that the second term on the r.h.s. vanishes for <p{t = 0) = 0. 
Backtransformation yields the following final equation for ip 

<p\t) = R{1 + Afi2) [3^ srnv? + (1 - 2/4) (cos <p - e 2 ^)] + e 2 ^ 2 (O) . (5) 
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It can be integrated numerically by separation of variable but the general class of the 
motion can already be seen by a stability analysis. By setting <p = <p = in the original 
equation (|l|) one sees that the fixpoint must fulfill the condition 

sin ip 

= tan w = a w 

COS if 

which is nothing but the initial definition for <po . To study the stability of the obtained 
fixpoint, one sets tp(t) = ipo + s(t) in Eq. ([!]) and neglects higher terms in e and e. The 
equation of motion for e reads 

9 

e = -— £ . 

H cos ipo 

This states that the fixpoint is marginally stable and the general motion is an oscillation 
around it of the form 

<p(t) = <p + Ae iut + Be- iut (6) 

where 

u:= ! 9 

y R cos tpo 

This oscillation can be seen in the simulation result shown in Fig. |l]b. It is also found for 
the motion of a particle on a plate M . Please note that the frequency of the oscillation 
in Eq. (||) does not depend on the external angular velocity, Q, as long as the full-sliding 
condition is fulfilled. 



3 Numerical Model 

In order to study the full dynamics of particles in a rotating drum, we use two-dimensional 
molecular dynamics (MD) simulations. Particles are approximated by spheres and only 



contact forces during collisions are considered jlljl . We use a linear spring- dashpot model 
for the force in normal direction, reading 

F n = ~k n h - y n h (7) 

and a viscous shear force in tangential direction, reading 

F s = -sign(7 g t> re i)min(7 s |w rel |, j u|F n |) . (8) 

Here h stands for the virtual overlap of the particle with the wall or with another particle, 
fi denotes the friction coefficient and v ie i the relative velocity of the two surfaces. The 
Coulomb friction law states that the magnitude of the shear force can at most have a 
value equal to the magnitude of the normal force multiplied by the friction coefficient. 
The model parameter k n is related to the material stiffness, j n to the normal restitution 
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coefficient, e n , and 7^ is chosen high enough to ensure that mostly the Coulomb friction 
law applies. 

Using MD simulations, the angle <p as function of time for a typical run is shown in 
Fig. |b using the parameters $7 = 1 Hz, e n = 0.9 and a particle-wall friction coefficient 
of [i w — 0.47 with the initial condition (p(0) — <p(0) — 0. The oscillatory motion given 
by the analytic result, Eq. @, can be seen as well. 

4 Results and Discussion 

In order to determine the maximum amplitude of this motion, Aip, one calculates the 
second turning point of the motion by setting tp = in Eq. (||) and using the initial 
conditions p(0) = and (fi(0) = 0. One obtains 

sin = 1 ~ 2 t& ( c 2m» a v _ CQS a^) , Aip > . (9) 

The sin-function on the l.h.s. can take only values from -1 to +1 and the expression in 
parentheses on the r.h.s. is always greater than zero which gives for small values of \i w 
a solution close to zero. The fraction on the r.h.s. changes sign at \x w = 1/V2 ~ 0.7071 
leading to a diverging solution around this value, e.g. for \i w = 0.71, Eq. @ does not have 
a solution. For a more quantitative analysis, the equation was integrated numerically 
using Mathematica and the obtained solution is shown as solid line in Fig. ||a as 
function of the friction coefficient fi w . Even though solutions of A<p > 90° are possible 
analytically from Eq. (||), they do not lead to an oscillating motion since for <p > 90° the 
particle will detach from the wall on its way downwards which invalidates the assumptions 
leading to Eq. (|l|). We added in Fig. ^a the data for Aip from the MD simulations as 
open circles (o). These were obtained by measuring the angle difference via the maxima 
and minima in plots like the one shown in Fig. [l]b. The values vary very little in time 
and the agreement with the analytic expression, taken from Eq. (^), is perfect, thus 
demonstrating that the simulation indeed describes the full-sliding regime when a viscous 
shear force is used together with the Coulomb threshold criterion, Eq. (||). For angles ip 
greater than 90 degrees, i.e. points lying above the dotted line in Fig. |2|a, the value of 
Aip was obtained by taking the maximum value during the first oscillation just before 
the free fall during the downwards motion sets in. 

The maximum amplitude is only reached for a high enough value of the rotation 
speed of the drum, f2, for the chosen initial conditions. For lower values of the system 
will oscillate with an intermediate amplitude since they are all marginally stable. If 
the initial conditions are chosen such that the initial amplitude is greater than Aip, the 
system will dissipate energy and finally reach a state with an amplitude lower than or 
equal to A<p. In the numerical simulations, the frequency of the oscillation shows no 



5 



180 - 



120 - 



1 1 1 
1 1 1 1 


1 1 

1 If 


analytic 




simulation o 

- 


f - 




(a) - 


i i i 


i i 



0.2 



0.4 



0.6 



35 
30 

<¥>> (°) 
25 



20 



15 



1 

• • 


1 1 

)Ste — — — e e — 


— 


— ; • & 

* y 


• 




* w 






; • 




(b) - 


if 


experiment ■ 
simulation 


■ - # — 

— e — 


• 

1 


i i 





10 

fl (Hz) 



15 



20 



Figure 2: (a) Comparison of the analytic expression for the maximum amplitude with 
the simulation and (b) frequency dependence of the average angle for the 
experiment (•) and the simulation using fi w = 0.47 (o). 



dependence on the rotation speed in agreement with the analytic result, Eq. (||), as 
long as the full-sliding condition is fulfilled. The frequency is 16.5 Hz for the simulation, 
extracted from Fig. [TJb, and 16.9 Hz for the analytic expression, Eq. (^). A similar system 
was investigated experimentally || by gluing two spheres together and slightly tilting the 
drum in order to force the particle to slide along one of the end caps. There, only a weak 
frequency dependence was found for values of £1 > 2.5 Hz giving a dominant oscillation 
frequency of 15.7 Hz, but for lower values, a stick-slip motion was clearly visible and 
led to lower values of the dominant frequency mode in the oscillation. Since no sticking 
was considered in the analytic model, it cannot capture this feature. The numerical 
model assumes a perfectly shaped drum whereas in the experiment a slight aspherity 
was observed and additional friction with one end cap was present thus amplifying the 
stick-slip region. Keeping these differences in mind the agreement is very satisfactory. 
By looking at the mean angle, (ip), of the particle as function of f2, a similar frequency 
dependence was found in the experiment. To a certain extent, this is also visible in the 
MD simulations and we show in Fig. ^3 the experimental points, taken from ||, as filled 
circles (•) and the numerical data by open circles (o). For O, < 7 Hz, the mean angle 
decreases with decreasing angular frequency of the drum and the general trend is well 
captured by the simulation. The large fluctuations in the experimental data can probably 
be attributed to the stick-slip motion which was enhanced by a slight aspherity of the 
drum. The best agreement with the experimental data in the high angular frequency 
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Figure 3: (a) Transition frequency to the centrifugal regime as function of friction 
coefficient and (b) angle of the particle as function of time for \x w = 0.665 
and f2 = 43 Hz (regime II). 



regime was obtained for fi w — 0.47. 

The solution of the analytic Eq. (|9|) and the corresponding Fig. |2|a state that the 
amplitude of the oscillation will be larger than 180° (= ir) for values of fx > 0.7 for a 
sufficiently high value of O. Physically that means that the particle will make overturns 
following the external sense of rotation and touches the wall all the time. This regime 
will be called the centrifugal regime and marked by the symbol III. It exists a well defined 
transition frequency which we will call £l c . The question now arises under what conditions 
this centrifugal regime exists and if it will be possible to make immediate overturns 
starting with a particle at rest. To address this question, we investigated numerically 
the transition frequency, J7 C , as function of friction coefficient, //, for a particle initially 
at rest. It can be viewed as a phase diagram and is shown in Fig. ||a. Three different 
regimes can be identified: (I) regardless of the angular frequency and the initial conditions 
the final state will oscillate below an angle of 90 degrees showing no overturns. (II) 
above a threshold of [i ~ 0.65, the centrifugal regime can be reached and the transition 
frequency decreases with increasing value of fi. Since the dynamics depends on the 
initial state, fluctuations in O c are visible. To illustrate this point, we show in Fig. |5^b a 
typical trajectory, particle angle (p as function of time, in this regime for the parameters 
fi = 0.665 and O = 43 Hz. After undergoing oscillations with different amplitudes 
initially, the particle reaches the centrifugal regime for t > 1.7 s. (Ill) for values of 
[i > 0.7 the particle will reach the centrifugal regime within the first external rotation. 
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The transition frequency is independent of the friction coefficient in this regime. 

The particle dynamics can be understood in more detail by looking at the motion in 
phase space, (ip, (p), see also ref. [12]. This is shown in Fig. |I]for different rotation speeds, 
f2, of the drum, ranging from 1 to 20 Hz and a particle starting initially at (ip, ip) = (0, 0). 
To the left, Fig. ||a, the motion in the variable ip is shown whereas to the right, Fig. |||b, 
99 is made dimensionless by dividing by f2. For small values of O, the particle initially 
sticks to the surface and rotates with the constant value ip = Q, leading to a horizontal 
line in phase space. For higher values of f2, a longer and longer transition period is 
visible and for the highest value, = 20 Hz, the particle never reaches the value tp = O 
but stays in full-sliding motion all the time. After the transition period, the particle 
oscillates around the value ip = arctan(^ w ) which shows up as circular motion around 
the point (<p,<p) = (arctan( / u W) ), 0) in the phase space plots. For our initial conditions, 
the maximum extend of the circular regime will be ip € [0°,90°] and (p 6 [— fi, fi] which 
is roughly given by the line corresponding to Q = 11 Hz in Fig. ||[ Please note that the 
centrifugal motion will correspond to a horizontal line at (p = SI. 

As discussed above, a simple regularized Coulomb friction law, see Eq. led to 
two different kinds of trajectories depending on the initial condition of the system. For 
initial amplitudes below A(p, the system will oscillate with this amplitude whereas for 
initial amplitudes above Aip, the system dissipates energy to reach a final oscillation 
with an amplitude lower than or equal to Aip. By considering different friction laws, e.g. 
choosing a different friction coefficient for the low and high velocity regime, up to four 



S 




different kinds of trajectories were found analytically [12 1 . The authors could also partly 
account for the stick-slip motion and achieved a closer agreement between theory and 
experiment. 

When a round particle, e.g. a sphere, is put in a rotating drum, it will start to 
rotate with the same sense of rotation as the drum and one might expect to be able 
to measure indirectly the coefficient of rolling friction by a similar method as described 
above. Unfortunately, the surface properties in our numerical model are too smooth and 
the angle of the particle, measured with respect to the vertical, will always give a value 
around zero, regardless of the angular velocity of the drum or the friction coefficient. 



5 Few Rotating Particles in a Two-Dimensional Drum 

If more than one round particle is placed in a two-dimensional rotating drum, the par- 
ticle rotations can be frustrated, i.e. suppressed, whenever particles have more than one 
contact point. The collective behaviour of all particles will then be very similar to that 
of a non-rotating particle as described in section |2] and the agreement becomes better 
with increasing particle number. 

To investigate the transition to the centrifugal motion in this case, a two-dimensional 
MD code is used which includes particle rotations. For the normal forces, the linear 
force with a linear damping was used and for the tangential force, the viscous friction 
law with a high enough parameter j s , to ensure that mostly the Coulomb friction law 
applies, was chosen. Two different friction coefficients will be used, namely: i) one for 
particle-particle collisions denoted by /i and ii) one for particle-wall collisions denoted by 
fj, w which will be varied independently. All particles have a diameter of d = 1 cm and 
the drum diameter was D = 10 cm. This gives a maximum of 28 particles that can form 
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Figure 6: (a) Average angle of 5 rotating particles for 17 = 50 Hz, /i = 0.75 and fj, w = 
0.35 and (b) critical particle number for the transition to the centrifugal 
regime as function of wall friction for three different particle-particle friction 
coefficients \i. The theoretical curve corresponds to the case fi = 0.00. 



the outer ring in the centrifugal regime and we limit our investigations to this value in 
order for all particles to feel the wall friction. 

For only a few number of particles, they will all form a line, shown in Fig. |5|a for 
N = 5, and the motion will essentially be one-dimensional. For more particles, where 
the exact number depends linearly on the drum to particle diameter ratio, D/d, more 
than one layer or ring of particles will form, shown in Fig. ||b for N = 15, and the motion 
is no longer one-dimensional but two-dimensional. Above a critical particle number, N c , 
which depends on D/d as well but also on the friction coefficient, fj,, all particles will 
form a single ring again but make full overturns jL3|. Such a situation is depicted in 
Fig. ||c for N = 20 and will be called the centrifugal regime in this case. 

First we look at the angle with respect to the vertical averaged over all particles, 
{ip), as a function of time. For a small number of particles, it shows a similar behaviour 
as the angle of a non-rotating particle which is shown in Fig. |6|a for the configuration 
from Fig. |5]a. Initially, the five particles form two layers but reorganize into one layer 
at t ~ 1.25 s which can be seen in Fig. |6|a by the change in amplitude. After this 
reorganization, oscillations around arctan(// w ), shown as dotted line, will occur as was 
the case for one non-rotating particle, see Fig. |l|b. 

To investigate the parameter regimes of the different motions depicted in Fig. [5], 
especially the transition to the centrifugal regime, we determine the minimal number of 
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particles, N c , which are necessary to reach it. For an inter-particle friction coefficient of 
/x = 0.0, all particles will form a single ring since piles cannot be supported. The particle 
configuration will fluctuate around a mean angle of (ip) = arctan(/i w ), see Fig. ^a. As 
discussed in the preceding section, a necessary condition for overturns is that the center 
of mass of the top-most particle reaches an angle of more than \. The angle 99 is 
measured from the vertical pointing downwards again, as shown in Fig. |l]a. Assuming 
that all particles are forming a single block and are touching their neighbours, the center 
of mass of the top-most particle of a configuration of A particles is at an angle of 
arctan(/i„,) + 75 (A — 4). This gives as necessary condition for overturns 

. , d , 1 . 7T 

arctanO^) + — (A - -) > - 
which when solved for the critical particle number, A c , gives 

A c > — (~ - arctan^)) + i . (10) 

This analytic expression for A c , derived for an inter-particle friction coefficient of 
fi = 0.0, is shown in Fig. ^b as function of the wall friction coefficient, fi w , as solid line. 
Also shown are the numerical results for fi = 0.25 (o) and fi = 0.50 (•) with a dotted 
line as least-square power law fit to the data points for the latter case to guide the eye. 
The differences with the analytic curve can be explained in the following way: If the 
inter-particle friction is sufficiently high, already two rotating particles can behave as 
a non-rotating particle in the two-dimensional drum. From the preceding section, one 
recalls that overturns are possible for one non-rotating particle if [i w > 0.65. Therefore, 
the critical particle number will decrease with increasing inter-particle friction for a high 
wall friction coefficient and we get lim^i A c = 2 for /j, w > 0.65. This is in perfect 
agreement with the numerical data where the data points for fi = 0.25 lie very close to 
the analytic curve for fi w > 0.45 and A c decreases with increasing /x giving a value of 
A c 5 for fj, = 0.5 and fi w = 0.65. 

For low wall friction coefficients, a larger number of particles is needed to reach the 
centrifugal regime since the value of arctan(/i„,) is small. If jx > 0, the particles can be 
arranged in more than one layer which gives a lower height for the top-most particle 
as compared to the one layer case discussed above. Consequently, A c will increase with 
increasing inter-particle friction for a low wall friction coefficient since the layer structure 
becomes more stable which is supported by the numerical data. 

6 Conclusions 

We investigated analytically the motion of a non-rotating particle in a rotating drum 
in the full-sliding limit. The angle of the particle's center of mass oscillates around a 
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mean value, given by the particle-wall Coulomb friction coefficient. The one-dimensional 
description breaks down for angles above 90°, corresponding to free flights of the particle. 
In order to study the full dynamics, including the transition from sticking to sliding 
motion for low velocities of the contact surfaces, we simulated the system using molecular 
dynamics. Full agreement with the analytic results in the full-sliding case was achieved 
and the general trend of the frequency dependence of the average particle angle observed 
in experiments could also be reproduced. The phase diagram of the transition to the 
centrifugal motion revealed three different regimes. 

We also found that the collective motion of a few rotating spheres in a two-dimensional 
rotating drum is, to a certain extend, similar to the motion of one non-rotating particle 
in a drum due to the frustrated particle rotations. We derived an analytic expression for 
the minimal number of particles to obtain centrifugal motion in the limit of negligible 
inter-particle friction. It compares well with the experimental findings and the change 
when inter-particle friction is added can be understood as well. 
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